Scar Intensity Statistics in the Position Representation 
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We obtain general predictions for the distribution of wave function intensities in position space 
on the periodic orbits of chaotic ballistic systems. The expressions depend on effective system size 
A'^, instability exponent A of the periodic orbit, and proximity to a focal point of the orbit. Limiting 
expressions are obtained that include the asymptotic probability distribution of rare high-intensity 
events and a perturbative formula valid in the limit of weak scarring. For finite system sizes, a single 
scaling variable A A'' describes deviations from the semiclassical N ^ oo limit. 
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I. INTRODUCTION 



Improving our understanding of wave function structure and transport in quantum systems with a non-integrable 
classical limit, and of their relationship with the dynamics of the corresponding classical system, has been a major focus 
of quantum chaos research. An early and visually dramatic finding along these lines was the presence of anomalously 
large fluctuations in wave function intensity near the short unstable periodic orbits of the corresponding classical 
^T^' dynamics. In 1983, McDonald observed individual wave functions (later called "scarred states") with a large 
' excess of intensity near a particular periodic orbit, as compared with the background Gaussian random fluctuations; 
■ the completeness condition then requires other wave functions (the "anti-scarred states") to have deflcient intensity 
near the same orbit. Since that time, scars have been recognized as the one of the few visually distinctive non-random 
features of quantum chaotic wave functions, as well as possibly the leading correction to the random wave conjecture 
made earlier by Berry and Voros 0, which states that quantum wave function statistics should converge to the 
behavior of random superpositions of plane waves in the semiclassical limit of short wavelength. 

The original theory of scarring was developed by Heller who placed a lower bound on scar fluctuations, as 
measured by overlaps of wave functions with Gaussian wave packets centered on the periodic orbits. Subsequent 
theoretical successes included work by Bogomolny 0|, who first computed an expression for scar intensity in the 
position representation (where scars were originally observed) , and by Berry Q , who developed a parallel analysis for 
the Wigner functions in phase space and predicted the structure of fringe patterns decorating each orbit. In these 
, pioneering investigations, the focus was mostly on average intensities for an energy window of nearby eigenstates, 
• rather than on properties of individual wave functions. Later, Agam and Fishman Iq] developed a semiclassical 
^n ' criterion for predicting scarring or anti-scarring in individual wave functions. In the semiclassical limit, the phase 
<^ : space region affected by the scar of any given periodic orbit tends to zero while the intensity enhancement in that 
^ ' region for a typical wave function remains finite; the scar effect is therefore fully consistent with the quantum ergodicity 
'rl^ , theorems of Shnirelman, Zelditch, and Cohn de Verdiere . 
C ■ Scars have been predicted and observed both experimentally and numerically in a wide variety of quantum systems 
J> \ with a chaotic classical limit (as well as in classical wave systems with a chaotic ray limit). Examples include 
r'"! ' applications as diverse as microwave cavities hydrogen atoms in a magnetic fiel d [p , electrons in a resonant 
rN \ tunneling diode with magnetic field J^J, Faraday surface waves vibrating soap films |l2l |. molecular vibrations 13], 
and acoustic radiation [l4]. 

In quantum or classical wave applications where dissipation effects are small enough for individual eigenstates to 
be resolved, it is of primary interest to study the scar effect for individual wave functions, rather than for averages 
over an energy or frequency window. In fact, the early empirical evidence for scarring as well as much of the empirical 
evidence since has been provided in the form of images of individual eigenstates. Such evidence, however, needs to be 
interpreted with care in light of the fact that random waves also exhibit regions of high intensity that to the eye look 
like scars but are in reality perfectly consistent with Gaussian random amplitude fiuctuations |l5j | . To understand the 
effect of scarring on physical processes in specific systems, and even to verify the presence of this phenomenon in the 
first place, it is necessary to have robust quantitative predictions concerning the statistical properties of individual 
wave functions, which can be compared with data as well as with the Gaussian random model. 

Antonsen et al. |l6j| noted that scar statistics of individual wave functions could be understood by combining 
known short-time dynamics associated with a given periodic orbit and random long-time recurrences whose statistical 
properties mirror those of random matrix theory (RMT). Similar ideas were used by Kaplan and Heller to obtain 
predictions for mean squared intensity on an orbit .17j , and also for the full distribution of wave function intensities on 
or away from the orbit [18| . These predictions were obtained in a Gaussian (Husimi) phase space basis, where classical- 
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quantum correspondence is most clearly manifest. On the other hand, wave function structure in position space is often 
of greater physical interest. Indeed, the original empirical evidence for scarring (Jj and much subsequent numerical and 
experimental investigation of scarring has focused on position-space properties, in contrast with the above-mentioned 
theoretical work. Wave function intensities in position space are directly accessible, for example, in microwave cavity 
experiments, and in the case of a diffusive billiard, comparisons have been made with predictions based on nonlinear 
sigma models 19] . Such explicit comparisons between observed deviations from RMT and theoretical predictions are 
lacking in the ballistic case. Furthermore, a quantitative understanding of position-space wave function statistics is 
essential for analyzing experiments on nanostructures such as quantum dots. For example, these statistics determine 
the effects of electron-electron interactions on conductance peak spacings in the Coulomb blockade regime 20] . 

This paper is organized as follows. The inverse participation ratio, or mean squared wave function intensity, on a 
periodic orbit is first obtained in Section^ followed by an analysis of the full wave function intensity distribution in 
Section IlIII In both cases, accurate and robust expressions for these statistical properties are obtained that depend 
only on the system size, instability exponent of the periodic orbit and proximity to a focal point. Several limiting 
expressions are obtained, including a result for the probability of rare high-intensity events, and a perturbative formula 
valid for moderate intensities when deviations from RMT are small. Finally, in Section Hvl we compare these results 
with data obtained in numerical simulations, and briefly investigate the behavior of wave function statistics beyond 
the semiclassical regime. 



II. INVERSE PARTICIPATION RATIO 



We begin by considering an unstable periodic orbit of period T in d dimensions, passing through periodic point 
(qotPo)- The Van-Vleck-Gutzwiller semiclassical propagator |2H includes a contribution from this periodic orbit to 
the return amplitude at time t = nT for position state qg: 
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Here S is the classical action associated with one traversal of the periodic orbit, and fi is the Maslov index, absorbing 

1/2 

changes in the sign of the amplitude. The amplitude ]det dp{qo, qt) / dqt\q[=qQ , where p is the initial momentum needed 
to travel from fixed go to qt in time is evaluated for the return trajectory qt = qo, and is the square root of the 
corresponding classical focusing factor. Of course, qo-, p, and qt are all d— component vectors in the general case. 
Omitted terms in Eq. are associated with other classical paths of length t beginning and ending at the same point 
qo; as we will see below the periodic path dominates the statistics of individual wave functions when the instability 
exponent of this periodic path is small. For sufficiently short times t or sufficiently small h, G^'^ {qo , qo , t) is a good 
approximation to the true quantum return amplitude G{qo,qo,t). 

For simplicity of presentation, and for easy comparison with numerical results in Section IIVI we now restrict 
ourselves to the case of discrete-time maps, which are commonly used as models for more general chaotic dynamical 
systems |22j . One-dimensional maps may be obtained by stroboscopical ly v iewing a one-dimensional system with 
time-dependent Hamiltonian (such as the kicked rotator or standard map |23), or alternatively by taking a Poincare 
surface of section for a two-dimensional time-independent Hamiltonian (such as a hard-wall billiard or smooth two- 
dimensional potential). For a one-dimensional map on a compact phase space of area 1, the quantum vector space is 
spanned by a finite basis of = l/27r?i independent position states \qi). Conventionally, we adopt a wave function 
normalization where the average of the discrete wave function intensities is set to unity: X]i3=o^ — the 

propagator G{qt,q,t) becomes a unitary N hy N matrix. With this normalization, 
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G^^ap(9o, 90, t) = { {^Nvn dq, ^^^^^ - ^ u ^ ^2) 

1 n = 

where n = t/T is an integer. Furthermore, for a, d = \ map, or equivalently for a d = 2 continuous-time dynamics, 
the monodromy matrix of the unstable periodic orbit is simply a 2 by 2 matrix, and may be written as 

dqt{q,p)ldq dqtiq,p)/dp \ _fa b \ / e"^ \ f d -b\ 
dptiq,p)/dq dptiq,p)/dp J - {c dj{ e""^ j{-c a J ^ 



where all partial derivatives are evaluated at the periodic point {q,p) — (goiPo), area preservation implies ad—bc= 1, 
and A > is the dimensionless positive instability exponent for one iteration of the periodic orbit. Then the inverse 
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of the focusing factor in Eq. ^ becomes 

{\dpiqo,qt)/dqtU,=go)~' = |5gt(9o,p)/5p|g.=go = 2|a6 sinh(AVr)| . (4) 

We are interested in the effect of the unstable periodic orbit on individual wave function behavior at or near qo, and 
specifically on the wave function intensities |'I'(go)P- A general discussion of the intensity distribution 7'(|5'(qo)P) is 
deferred until Section IlIII here wc focus first on the mean squared intensity, also known as the inverse participation 
ratio (IPR), 



i(9o) = l*(go)l^ 



(5) 



where the average is performed over eigenstates of the system, and the mean intensity |^'((j'o)P is normalized to 
unity as above. For a general Hamiltonian system, the average must be restricted to a classically narrow energy 
window; in the case of a discrete-time map, however, we are free to average over all values of the (periodic) quasi- 
energy. T^q^) is the inverse fraction of wave functions in the energy window that have significant intensity at q^; its 
possible values range from 1 in the case of complete ergodicity (where every wave function has equal intensity at go) 
to N in the case of complete localization (where the wave functions are delta functions in position space). 

As has previously been discussed 24] , for a non-degenerate system the IPR is proportional to the long-time return 
probability for initial state \qo}, which for the chaotic case is in turn proportional to a sum of return probabilities due 
to the periodic orbit. 



^{qo) = / dt|Gmap(go,9o,OI' 



F 



Sr=-oo|gE.ap(go,go,»r)l^ 
|Gg;ap(go,9o,0)P 



(6) 



Here 0^°^ is the contribution to the propagator from a single dominant periodic orbit, as given in Eq. J^l, and A/" is a 
normalization constant. The multiplicative factor F encodes information about all recurrences not included in GP;°p; 
semiclassically these recurrences are associated with closed classical trajectories other than the original periodic orbit. 
Treating these other contributions as random and uncorrelated, which is physically justifiable for chaotic dynamics at 
long times, leads to the random wave result F = 3 (for real wave functions, in the presence of time reversal symmetry) 
or = 2 (for complex wave functions, in the absence of time reversal symmetry). 
Combining the results of Eqs. ©, and Q), we obtain the prediction 
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For instability exponent A 3> 1, only one iteration of the orbit contributes, and we obtain Ipiod (go) = F 1 
converging eventually to the random wave result XRMT(go) = F. In the more interesting limit A ^ 1, we obtain 
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where 1 <C M ^ A ^ and 7 « 0.577 is Euler's constant. Then we find 



2'prcd(go) — F 
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in the weakly unstable regime A ^ 1. 

From Eq. Q or Eq. (j^, we see that scar strength, as measured by the inverse participation ratio, depends on three 
parameters only: the semiclassical parameter N (or ratio of system size to wavelength), the instability exponent A 
of the periodic orbit, and the phase space orientation parameter |a6|, which is directly related to the angle between 
the p— axis and the unstable manifold of the orbit at the point (q'Oj^'o)- In a two-dimensional Hamiltonian system, 
such as a billiard or smooth potential, the stable manifold rotates and the \ah\ parameter changes as one moves along 
the orbit, approaching zero at the focal points. Thus, all possible values of \ah\ are relevant for describing the wave 
function behavior near a generic orbit, while A is a constant parameter for a fixed orbit. For small values of |a6|, 
which correspond to the strongest scarring, \ah\ is directly proportional to the distance from the nearest focal point 
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on the orbit. In a specific discrete-time map, \ab\ has a single value for each orbit of period one, since each such orbit 
has only one periodic point. However, the same correspondence between the \ab\ parameter and distance to a focal 
point applies also in the case of maps, if we consider a family of maps arising from Poincare sections intersecting a 
given orbit at various distances from a focal point. 

We need to understand the range of validity for Eq. lO or Eq. (jSJ. Classical-quantum correspondence G{qo, qo,t) ~ 
G^'^{qo,qo,t) fails when qo is within a wavelength of a focal point, as the classically singular probability density is 
smoothed out by the uncertainty principle. Thus, our semiclassical derivation fails for \ab\ < N~^, and the correct 
scaling behavior for this regime must be obtained by replacing |a6| with a number of order unity. Thus, the inverse 
participation ratio X(go) remains finite at 0(A~^) for A ^ 1 as |a6|A^ 0. We note that the IPR scales as A^^ in the 
basis of phase-space Gaussians optimally oriented with respect to the stable and unstable manifolds of the periodic 
orbit |l7j : it is not surprising to obtain the same scaling behavior in position space near the focal points of the orbit, 
since a position state in this case may be thought of simply as a limiting case of a family of optimal Gaussians. 

A second limitation on the validity of Eq. lO or Eq. Q is that the Lyapunov time T/A needed to escape from the 
vicinity of the orbit must be much shorter than the Heisenberg time N, at which quantum mechanical exploration 
ceases and the quantum dynamics becomes quasi-periodic. For AA'^ < T, the infinite sum in Eq. IJJ) must be cut 
off at n ~ N/T to avoid counting classical recurrences that occur after the Heisenberg time and have no quantum 
analogue. This has the effect of replacing AA^ in Eq. (j^ with a constant of order T when AA^ < T, so that the inverse 
participation ratio again remains finite: T{qa) ~ \ab\^^T^^ as A 0. The behavior of the IPR for AA^ < T will 
be discussed further in Section Hvl For A and \ab\ simultaneously, the IPR must obey an upper bound 
^{qo) < which is saturated only if the position state |go) is itself an eigenstate of the dynamics. 

In conclusion, for a given semiclassical parameter A^, we require \ab\ > and A > TN'^ for our semiclassical 
expressions to be valid; smaller values of A or \ab\ do not lead to parametrically stronger scarring. At the same time, 
as seen from Eq. (Q, we need simultaneously A <C 1 and \ab\ <C 1 to obtain a large enhancement of the IPR above 
the RMT prediction. We also note that in any specific dynamical system, multiple trajectories (periodic or closed) 
contribute to I{qo) at order A^"^, and a single orbit can only be expected to dominate the wave function statistics at 
qo when \ab\ <Si 1 and A ^ 1, i.e. when qo is close to a focal point of a weakly unstable periodic orbit. 



III. SCAR INTENSITY DISTRIBUTION 



Although the IPR provides a useful one- number measure of the degree of wave function localization at qo, it is only 
the lowest non-trivial moment of the intensity distribution, and as such provides limited information about the full 
distribution. In particular, T{qo) > Irmt = F implies an increased intensity variance and clearly suggests longer tails 
of the intensity distribution, i.e. an enhanced probability of finding wave functions with intensity / = |^(go)P much 
larger or much smaller than the average intensity. However, to obtain quantitative predictions about the probability 
of such rare events, we must go beyond the IPR to examine the full intensity distribution 

As in Section ^ we assume a clear separation of scales exists between known short-time recurrences associated 
with a particular short periodic orbit and new long-time nonlinear recurrences which semiclassically may be associated 
with an exponentially large number of homoclinic paths starting and ending near qo. The statistical properties of 
these new recurrences at long times may be assumed to be consistent with RMT, i.e. independent of the fact that qo 
happens to lie on a short periodic orbit 17]. The full return amplitude at long times then becomes a convolution of 
the known short-time dynamics and the RMT- like random long-time recurrences [T^ : 

G„,ap(go, 90, <) « E GE.^plgo, 90, r)G™'i(Qo, 9o, t - r) . (10) 

T 

Fourier transforming, we find that the local density of states at qo is given by the product of a smooth envelope 
associated with the periodic orbit and a random fluctuating part: 

l*a(g)|' SiE - E^) = G^ap(go, qo, E) = GP,:p(go, go, i?) ^ \ro.?5{E - E^) , (11) 

a Q 

where the Tq, are distributed as (real or complex) independent Gaussian random variables with mean and variance 
1, and GP;°p((j'o, 9o, E) is the Fourier transform of Eq. 101. The individual wave function intensities are given by 

= |*a(go)|' = G^:p(go,go,i?a)kaP = S{E^)\r^W (12) 
where we have given the name S{E) to the smooth part of the local density of states at qo- Explicitly, this smooth 
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envelope coming from the short periodic orbit takes the form 

S{E)^l + 2^^± , ' cosf '^^'^ + ^^ ^^-^V (13) 

The fuU intensity distribution at go is given by an energy average. Since the level density is constant for a map (and 
constant within the energy window of interest for a continuous-time Hamiltonian system), and the random factors 
jrctp are uncorrelated with the smooth envelope S{E), each moment of the intensity distribution may be obtained as 
a product of a factor depending on properties of the orbit and a universal factor associated with Gaussian random 
fluctuations: 



1 l-Ei+2TTh/T \ 

dE {S{E)r\ ■ (|r„P«) . (14) 



27rh/T 



El 



Note that the smooth envelope is_periodic i n T, so we need only to average over the energy interval Ei < E < 
El + 2'Kh/T. In particular, using P = ^{qo), Ifa]^ — F, and Eq. l(T^ . we recover the result of Eq. Q for the inverse 
participation ratio. 

In the absence of time-reversal symmetry, is a complex Gaussian and the random factor in Eq. (|12|) is exponen- 
tially distributed, 

^(Kn = e~l^°l% (15) 
while in the presence of time-reversal symmetry Tq, is real and its square follows a Porter-Thomas distribution, 

PTRsdraH = ^=4=^e-l'-°l'/2 . (16) 



For definiteness, we focus on the generic situation of no time-reversal symmetry. At a fixed energy E, the distribution 
of / = S{E)\ra\'^ is V{I) ~ 'slE)^^^^^^^^ ^ combined intensity distribution over all energies is 

V{I) = ^ / . (17) 

2'Kh/T J S[E) 

Several limits are of particular interest. If the scarring is relatively weak, i.e. S{E) — 1 <C 1 in Eq. H13|) . we may 
expand Eq. H17|l to second order in S{E) — 1 to obtain 

V{I) = (l + e - 2e/ + |/2 + 0(e2)) , (18) 

in complete analogy with perturbative results for disordered systems, in the limit of weak disorder l^^. Here 

. 1 r+''^'/^z..c..^^2_^M ,_I(<?o)~2 



2nn/T 



dE {S{E)Y = ^ _ 1 = , (19) 
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and for e — > we recover the Porter-Thomas distribution predicted by random matrix theory: 'Prmt(^) = . In 
the perturbative regime, deviations of the intensity distribution from the RMT prediction are proportional to the 
deviation of the IPR from its RMT value of 2. We note that the validity of Eq. H18|l requires / <C 1 in addition to 
/e <C 1, so perturbation theory does not apply in the tail of the intensity distribution. 

To obtain a simple expression for the probability of rare events, i.e. intensities much higher than the average 
intensity / = 1, we need to follow a different approach. For large /, the integral in Eq. H17|l may be evaluated using 
the saddle point (stationary phase) method, with the dominant contribution coming from the energy at which the 
smooth envelope S{E) is peaked. Physically, this means that very large wave function intensities will almost always 
be at the "scar energies" Eq where interference from successive iterations of the orbit add fully constructively, and 
never at the "antiscar energies" where this interference is destructive. We then obtain 



where Eq is the energy at which the smooth envelope S{E) of Eq. (|13|l is peaked, and S{Eq) and S"{Eo) are the value 
and the second derivative of S{E) at the peak. 
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For small A, we may rewrite the sum in Eq. H13() as an integral, which may be evaluated in terms of hypergeometric 
functions. We then find that the smooth local density of states is peaked at 

^ -S + finh/2 + f]Xh 



where the dimensionless constant /3 w 0.2133 is the solution of dx xsm((3x — 7r/4)/\/sinhx = 0. The height of the 
peak is given in this limit by 



while the second derivative is 



where 7 = dx cos{(3x — 7r/4)/-\/sinha; sa 3.059 and S = dx x"^ cos(/3x — 7r/4)/-\/sinhx « 16.14 are dimensionless 
constants. Combining these results, the tail of the intensity distribution for small A finally takes the form 



where A = ^J2/N\ah\ is a geometrical factor that becomes 0(1) near a focal point of the periodic orbit. Comparing 
with 'Prmt(-^) = e~^, we note that the tail of the intensity distribution remains exponential in the presence of scarring, 
but the exponent now depends on the instability exponent A, leading to a greatly enhanced probability of finding 
very high wave function intensities on weakly unstable orbits. The exponentially small probability of rare events in 
ballistic chaotic systems, as given by Eq. ()24|l . contrasts with the log-normal tail predicted and observed for diffusive 
two-dimensional systems in the metallic regime |25l |. 

In the opposite hmit / — > 0, Eq. l(T7ll implies V{I) — > > 1/ S — \ = V]xmt{I)', this enhanced probability of very 
small intensities is known as the antiscar effect. 



IV. NUMERICAL MODEL AND RESULTS 

For the purpose of testing the general results obtained in Sections ^ and IIIII we consider a specific family of 
discrete-time kicked maps defined classically on a toroidal phase space (g,_p) G [0, 1) x [0, 1). The time evolution for 
one step ((j„,_p„) ^ (q„+i,p„+i) is given by 

<ln+l/2 = qn+T[{pn) mod 1 

Pn+1 = p„ - V'{qn+i/2) mod 1 (25) 

Qn+l qn+1/2 +T2{pn+l) mod 1 

with 

Ti{p) - ^(p-po)^ + ^1 [2sin27r(p-po) -sin47r(p-po)] 

V{q) = -^Kg{q~qof +D2Kg[2sm2TT{q-qo)-sm4TT{q-qo)] (26) 

T^ip) = ^Kp{p-pof +D-iKp[2sm2TT{p-po)-sm4TT{p-po)] . 

Physically, such a map may be thought of arising from a time-periodic Hamiltonian where a particle experiences 
free evolution under the influence of kinetic term Ti{p) for the first half of a period, followed by a sudden kick of 
strength V{q) and then by additional free evolution with kinetic term T2{p) for the second half of the period. The 
corresponding quantum evolution for one step is given by 



(j ^ g-^T2{p)/h^-lV(q)/h^-lT,{p)/h ^ 



(27) 
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q and Kp (so that both V'{q) and T2{p) are continuous functions on the torus) would correspond 



where h = 1/2ttN, and N is the dimension of the Hilbert space. If we were interested in the spectral behavior only, 
we could of course view the same dynamics stroboscopically right before or right after each kick, effectively combining 
the two free evolution sub-steps into one governed by kinetic term Ti +T2. In that case, choosing = Po = and 
integer values for 

to a perturbed cat map, a system well studied in the literature (26j 

In our case, we choose non- integer values ot Kq and Kp, as well as generic values for qo andpo, resulting in diffraction 
at (7 = and p = 0. Nonzero qq and po also cause a breaking of time-reversal symmetry. By construction, the above 
map is fully chaotic, for positive Kq and Kp + 1, and has a periodic orbit of period T = 1 at {q,p) = {qo,po). The 
instability exponent A and unstable manifold orientation parameter \ab\ of this orbit are implicitly determined by the 
parameters Kq and Kp via the relations 



2coshA = 2 
2|a6|sinhA = 1 



Kq{l 



K„ 



Kp{l+Kq 



(28) 



Any desirable values of A and \ab\ may be obtained by selecting appropriate Kq and Kp in accordance with Eq. I|28|l . 
The parameters Di, D2, and have no effect on the unstable orbit at ((jOiPo) or on its monodromy matrix (Eq. Q), 
and within certain bounds these additional parameters do not affect the fully chaotic nature of the dynamics 26] . Such 
additional parameters may be used to construct a large ensemble of distinct classical systems with identical linearized 
behavior near the dominant periodic orbit. In our numerical investigations, however, we did not see significant effects 
in wave function statistics due to variation of D2, and -D3, and for simplicity the data shown below was collected 
using Di — D2 = -D3 — 0. For a given pair of parameters Kq, Kp (equivalently. A, \ab\), we do generate an ensemble 
of similar systems by varying q^ and Pq, as well as varying boundary conditions in both the q and p directions. 
Specifically, the boundary condition variation is accomplished by taking parameters ei and £2 uniformly distributed 
between and 1, and requiring '^{q -I- 1) = '^{q) e~*^^^i and '^{p -I- 1) = ^'(p) e*^'^^^; this is equivalent to the choice 
qj = {i + €2)/N, < j < N - 1 for the position-space basis and pe ^ {£ + ei)/N, < £ < N - I for the 
momentum-space basis. 
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FIG. 1: The average inverse participation ratio T{qo) on a periodic orbit is shown as a function of instability exponent A of the 
periodic orbit, for the system defined by Eqs. 12511 and 12611 with system size TV = 1024. The upper and lower curves correspond 
to \ab\ — 0.005, 0.025, respectively. The scar theory prediction of Eq. Q with _F = 2 is indicated by dashed lines. All quantities 
shown in this and subsequent figures are dimensionless. 



In Fig.^ we examine the IPR as a function of the instability exponent A of the periodic orbit, for two values of \ab\. 
Good agreement is observed with the simple asymptotic expression given by Eq. @ for a wide range of parameters, 
as the size of the observed and predicted fluctuations varies by an order of magnitude. We recall that the validity of 
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Eq. requires <C A <C 1, since the orbit period T = 1 in the present case. Fig. 13 shows what happens to the 
IPR as we leave the semiclassical hmit ^ A. We observe that to a very good approximation, deviations from 
the semiclassical prediction are well described by a simple scaling variable AA^: 

X(go) = /(AiV)2'prcd(go), (29) 

where 2prcd(9o) is the asymptotic prediction of Eq. f{x) ^ 1 for a; ^ 1, and f{x) ~ x for a; <C 1 (as indicated by 
the dotted line in Fig.EJ. The dimensionless parameter XN may equivalently be expressed as the ratio Th/T\, where 
Th is the Heisenberg time at which individual eigenstates are resolved and = 1/A is a Lyapunov time associated 
with classical decay away from the periodic orbit. It is noteworthy that the crossover between the semiclassical 
behavior at large XN and the saturated behavior at small AA^ occurs around XN « 40 for all parameter values 
considered. This implies, for example, that in applications to quantum dots in the Coulomb blockade regime, where 
30 < < 70 in some typical experiments |27l |. asymptotic N ^ oo expressions are inadequate for describing 
wave function intensity statistics in "generic" chaotic mean-field potentials (A ~ 1), and finite-7i expressions must 
be used instead. Accurate modeling of wave function statistical properties in such systems is necessary for proper 
understanding and interpretation of conductance peak spacing experiments |2Qj. 
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FIG. 2: f{\N) (Eq. lEHJ), the ratio of the observed inverse participation ratio to the asymptotic prediction of Eq. Q, is shown 
as a function of scaling parameter AA'^, for three different system sizes. Open and filled symbols represent \ab\ = 0.025 and 
0.05, respectively. The dotted line indicates the asymptotic behavior f{XN) ~ 0.035 AA^, for AA'^ <^ 1. 



We now turn to the full distribution of wave function intensities. A typical probability distribution is shown in 
Fig. 13 and the scar theory prediction of Eq. (|17|l clearly describes the data much better than does the random matrix 
theory result. We note that the distribution of rare high-intensity events is in adequate quantitative agreement with 
the simple "tail" formula of Eq. H24|l . although the latter result is strictly valid only in the A — > limit. In Fig. ^ 
we see that the scar theory prediction also works well in the region of small to moderate intensities; in particular it 
correctly predicts the antiscar enhancement for / ^ 1. The parameters used here are evidently outside the range of 
validity of the perturbative expression (Eq. For the values of TV and \ab\ considered in Fig. 01 the perturbative 

result becomes applicable in the main body of the distribution for A > 1.5 (not shown), where the system is sufficiently 
unstable that deviations from RMT fall below the 20% level. 
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FIG. 3: The wave function intensity distribution for A*" = 512, A — 0.9, and \ab\ = 0.005 is compared with the scar theory 
prediction of Eq. II17|I (solid line) and with the simpler expression of Eq. 1)24^ (dashed line), which is valid for large intensities. 
The random matrix result ■Prmt(/) = is shown by a dotted line for comparison. 




FIG. 4: The wave function intensity distribution ?-'(/) is scaled by the random matrix theory prediction 7'rmt(/), and the 
result is compared with the full prediction of Eq. 11711 and the perturbative limit of Eq. I18II . All parameters are the same as 
in the previous figure. 
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V. SUMMARY 

We have studied the distribution of position-space wave function intensities on an unstable periodic orbit in a 
ballistic chaotic system. In the two-dimensional case (eqiiivalently. in a one-dimensional discrete-time map), the full 
distribution of intensities is given by three parameters: (i) the system size in wavelength units, or alternatively the 
dimensionless conductance, inverse effective U, or Heisenberg time in units of a one-bounce time; (ii) the instability 
exponent of the periodic orbit; and (iii) orientation of the unstable manifold or alternatively the proximity to a focal 
point. Generalization to higher dimensions is straightforward. 

The tail of the probability distribution is exponential, with exponent proportional to a simple function of the 
above three parameters. This contrasts with the log-normal distribution of rare events in disordered two-dimensional 
systems. When the scar effect is weak, the behavior for moderate intensities is given by a perturbative expression in 
complete analogy with the disordered case. 

Simple semiclassical expressions begin to break down when the classical Lyapunov time associated with decay 
away from the periodic orbit reaches 2.5 — 3% of the Heisenberg time. A better quantitative understanding of these 
saturation effects is needed for reliable predictions of wave function statistics in finite-size ballistic systems. 
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